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O ' 

We present a theoretical study of the emission from a superluminal polarization current whose 
distribution pattern rotates (with an angular frequency u>) and oscillates (with a frequency Q) at the 
(— i , same time, and which comprises both poloidal and toroidal components. This type of polarization 

. current is found in recent practical machines designed to investigate superluminal emission. We find 

that the superluminal motion of the distribution pattern of the emitting current generates localized 
electromagnetic waves that do not decay spherically, i.e. that do not have an intensity diminishing 
£f>l . like Rp~ 2 with the distance Rp from their source. The nonspherical decay of the focused wave 

packets that are emitted by the polarization currents does not contravene conservation of energy: 
the constructive interference of the constituent waves of such propagating caustics takes place within 
different solid angles on spheres of different radii (-Rp) centred on the source. For a polarization 
current whose longitudinal distribution (over an azimuthal interval of length 2n) consists of m cycles 
of a sinusoidal wave train, the nonspherically decaying part of the emitted radiation contains the 
frequencies £l±muj; i.e. it contains only the frequencies involved in the creation and implementation 
of the source. This is in contrast to recent studies of the spherically decaying emission, which was 
shown to contain much higher frequencies. The polarization of the emitted radiation is found to be 
. ^ , linear for most configurations of the source. 

CZ2 
>v 

Oh' 



Or 
O 



> 

m 
oo 



o 



Oh 



- 1—1 

X 

5-H 



PACS numbers: 



I. INTRODUCTION 



The electromagnetic field of a moving charged particle whose speed exceeds the speed of light in vacuo is the subject 
of several papers written by Sommerfeld in 1904 and 1905 [l], papers which were regarded to have been superseded by 
special relativity soon after their publication. One reason for abandoning further investigations of the work initiated 
by Sommerfeld, of course, was that any known particle that had a charge also had a rest mass and so was barred 
from moving faster than light by the requirements of special relativity. But there was an additional reason. Not even 
a massless particle can move faster than light in vacuo if it is charged; for if it does it would give rise to an infinitely 
strong electromagnetic field on the envelope of the wave fronts that emanate from it. 

It was not until the appearance of the works of Ginzburg and his coworkers QEJQ that it was realized that, though 
no superluminal source of the electomagnetic field can be point-like, there are no physical principles disallowing faster- 
than-light sources that are extended. The coordinated motion of aggregates of subluminally moving charged particles 
of opposite sign can give rise to macroscopic polarization currents whose distribution patterns move superluminally. 
The electromagnetic field that is generated by such extended sources, on the other hand, may be built up by the 
superposition of the fields of their constituent volume elements, elements which individually act as the superluminally 
moving point sources considered by Sommerfeld [H IS IS @ • 

The purpose of the present paper is to examine the radiation field of a particular class of such volume-distributed 
sources: polarization currents whose distribution patterns have the time dependence of a travelling wave with a 
centripetally-accelerated superluminal motion. A motivation for this work is the recent design, construction and 
testing of experimental machines with this characteristic [jj 0, 0, 0, 0, 0] . 

The Green's function for the problem we analyze is the familiar field of a uniformly rotating point source. This is the 
Lienard-Wiechert field that is encountered in the analysis of synchrotron radiation, except that here we do not restrict 
the speed of the source to the subluminal regime. This uniformly rotating point source is used as a basic volume 
element of the extended superluminal source, which is then treated in detail by superposing the Lienard-Wiechert 
fields of its constituent elements. We find that fundamentally new radiation processes come into play as a result of 
lifting the restriction to the subluminal regime, processes that have no counterparts in either the synchrotron or the 
Cerenkov effects. 

This paper is organized as follows. Section II presents an introductory description of the principles behind our 
calculation of the emission. It includes a definition of frequencies (Table I) and polarization (Section HE) of the 
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source, which are later used to discuss the spectral characteristics and polarization of the emitted radiation. A 
detailed mathematical treatment is given in Section III, with the expression describing the nonspherically decaying 
component of the radiation [Eq. (57)] being derived and discussed in Section HID. The nonspherically and spherically 
decaying components of the emission are compared in Section IV, and a summary is given in Section V. 



II. PREAMBLE: FORMULATION OF THE PROBLEM 



The physical principles underlying the (mostly unexpected) consequences of the processes described in Section I are 
more transparent if we begin with a descriptive account of the building blocks of the (unavoidably lengthy) analysis 
that follows in Section III. Consider a localized charge q (a charge with linear dimensions much smaller than the 
typical radiation wavelengths considered) which moves on a circle of radius r with the constant angular velocity we z , 
i.e. whose path x(f) is given, in terms of the cylindrical polar coordinates (r, tp, z), by 

r = const., z = const., tp = tp + Lut, (1) 

where e 2 is the basis vector associated with z, and tp the initial value of tp. 

Having defined the path of what will constitute the basic element of our volume source, we shall use the remainder 
of this section to describe (A) the associated Lienard-Wiechert field, (B) the bifurcation surface which divides the 
volume of the source into parts with diffcrcing influences on the field, (C) the Hadamard regularization technique for 
dealing with the singularities which are inevitably encountered in the fields from a superluminal source, (D) the focal 
regions in which the nonspherically decaying part of the emission is detectable and (E) the relationship between the 
polarization of the source and of the emission. 



A. The Lienard-Wiechert fields 



The Lienard-Wiechert electric and magnetic fields that arise from a superluminally moving charge q with the 
trajectory described in Eq. (1) are given by 



E(x P ,tp) = qJ2 

tret 



(1 - |x| 2 /c 2 )(n — x/c) nx{(n — x/c)xx} 



|1 -h-±/c\ 3 R 2 (t) c 2 |l-n-x/c| 3 i?(t) 



(2) 



and B = nxE. Here, the coordinates (xp,tp) = (rp, tpp, Zp,tp) mark the space-time of observation points, x = 
dx/dt, c is the speed of light in vacuo, and R and n are the magnitude 

R(t) = [{z P - z) 2 + r P 2 + r 2 - 2r P rcos(tp P - tp - ujt)}? (3) 

and the direction n = R/i? of the vector R(t) = xp — x(t). The summation extends over all values of the retarded 
time, i.e. all solutions t TCt < tp of h{t) =t + R(t)/c = tp. This summation and the absolute-value signs in the factor 
1 — fi • x/c | , which stems from the evaluation of the Dirac delta function S(t — tp + R/c) in the classical expression for 
the retarded potential [Eq. (10) below], are omitted in most textbook derivations of Lienard-Wiechert fields because 
in the subluminal regime the retarded time is a single- valued function of the observation time and 1 — n • x/c is 
everywhere positive [lq . 

For a given point source (r, tp, z) with rto > c and various positions (rp,tpp,Zp) of the observation point, the 
dependence tp = h(t) of the reception time tp on the emission time t can have one of the generic forms shown in Fig. 
1. As can be seen from curve (a) of this figure, there are values t± of the retarded time at which 

1 - fi • x/c = h(t) = 1 - r P {ruj/c) sin{tp P - tp + uot)/R{t) = 0, (4) 

i.e. at which the source approaches the observer with the speed of light along the radiation direction fi. The set of 
observation points at which the field receives contributions from the values t± of the retarded time are those located 
on the envelope of the wave fronts emanating from the moving source in question (Fig. 2). In the vicinity of the 
extrema of curve (a) in Fig. 1, the transcendental equation tp — h(t) for the retarded time reduces to 

tp = t P± + \h{t ± )(t-t ± ) 2 + -.., (5) 

where tp± = h(t±) are the values of tp at which the waves emitted at t — t± arrive, and constructively interfere, at 
the envelope of the wave fronts (see Appendix C of 0). 



FIG. 1: The relationship between the observation time tp and the emission time t for an observation point that lies (a) inside 
or on, (b) on the cusp of, and (c) outside the envelope of the wave fronts or the bifurcation surface shown in Figs. 2 and 3. This 
relationship is given by tp = t + R(t)/c = h(t;r,<p, z;rp,ipp, Zp), an equation that applies to the envelope when the position 
(r, tp, z) of the source point is fixed and to the bifurcation surface when the location (rp, ipp, zp) of the observer is fixed. The 
maxima and minima of curve (a), at which dR/dt = — c, occur on the sheets (j>+ an d 4>- of the envelope or the bifurcation 
surface, respectively (Figs. 2 and 3). The inflection points of curve (b), at which d 2 R/dt 2 = 0, occur on the cusp curve of the 
envelope or the bifurcation surface (Fig. 4). 

Were it to exist, a superluminally rotating point source would therefore generate an infinitely large field on the 
envelope of the wave fronts that would emanate from it: the Lienard-Wiechert fields diverge at the extrema of curve 
(a) in Fig. 1, where the factor 1 — n • x/c in the denominator of Eq. (2) vanishes. However, superluminal sources 
are necessarily extended and so it is only a superposition of the Lienard-Wiechert fields of their individual volume 
elements that is physically meaningful. The relevant quantity is the integral of the above field over the volume of 
(r, tp, z) space that the moving localized source occupies in its rest frame. 

B. The bifurcation surface of an observation point 

Considering the field of a given point source, we have so far kept (r, tp, z) fixed and have varied (rp,<pp, zp). If 
we now keep the observation point fixed and allow the coordinates (r, tp, z) of the moving source point to sweep the 
rest-frame volume of an extended source, then the curves in Fig. 1 would represent the forms the relationship tp = h(t) 
assumes in different regions of the (r,(p,z) space. For any given (rp, tpp, zp), there is a set of source elements of a 
volume source which approach the observer along the radiation direction n with the speed of light at the retarded 
time, i.e. for which h — 1 — n • x/c is zero. The locus of this set of source points in the (r, tp, z) space is given by the 
intersection of the volume of the source with the two-sheeted surface tp±{r, tp, z; rp, tpp, zp) — tp = 0, a surface that 
we shall refer to as the bifurcation surface of the observation point P (Fig. 3). For all source points in the vicinity of 
this locus, the relationship between the retarded and the observation times has the form assumed by curve (a) of Fig. 
1 in the neighbourhood of its extrema, i.e. is that given in Eq. (5). 

The two sheets (±) of the bifurcation surface tp±(r, tp, z; rp, tpp, zp) — tp = meet tangentially along a cusp curve 
(Figs. 3 and 4). The source points on this cusp curve approach the observer not only with the wave speed (with 
dR/dt = — c) but also with zero acceleration (d 2 R/dt 2 = 0) along the radiation direction. For such source points, 
the relationship between the retarded and the observation times has the form assumed by curve (b) of Fig. 1 in the 
neighbourhood of its inflection point. 

The retarded times from which the Lienard-Wichert field (2) receives singular contributions during any given period, 
therefore, are t ~ t + ± [2(t P H — t P+ )/h(t + )]i for t P <t P+ and t ~ t_ ± [2(t P - t P -)/h(t-)]i for t P >t P _ [see Eq. 
(5)]. Close to these retarded times, the factor 1 — n • x/c has the absolute value 

|1 - n • x/c| = \h(t)\ ~ [2h(t ± )(t P - t P ±)]v (6) 
according to Eq. (5). The functions tp± = t± + R(t±)/c in Eq. (6) depend linearly on the coordinate tp of the source 
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FIG. 2: (a) Envelope of the spherical wave fronts emanating from a source point S which moves with a constant angular 
velocity uj on a circle of radius r — 2.5c/lo (r = ruj/c = 2.5). The circles in broken lines designate the orbit of S and the light 
cylinder rp — c/u (rp — 1). The curves to which the emitted wave fronts are tangent are the cross sections of the two sheets 
<f>± of the envelope with the plane of source's orbit, (b) Three-dimensional view of the light cylinder and the envelope of wave 
fronts for the same source point S. The tube-like surface constituting the envelope is symmetric with respect to the plane of 
the orbit. The cusp along which the two sheets of this envelope meet touches, and is tangential to, the light cylinder at a 
point on the plane of the source's orbit and spirals around the rotation axis out into the radiation zone. It approaches the cone 
8p — arcsin(l/f) as Rp tends to infinity (Rp, 9p and tpp are the spherical coordinates of the observation point P). 

point. This can be seen from Eq. (4) without solving for t±: since t in Eq. (4) appears in only the combination 
ipp — if — u)t, the solutions t± of this equation are given by expressions of the form t± — (ip± — tp)/u in which tp± are 
functions of (r, z) only. Given that the functions R(t±) also depend on t through ipp — (p — cut and so are independent 
of tp, it follows from tp± =t± + R(t±)/c that tp± —tp has the functional form (tp± ~ tp)/uj in which tp± depend on 
r and z only. 

According to Eqs. (2) and (6), therefore, the Lienard-Wiechert fields diverge like \tp± — tp\~z for those source 
elements close to the bifurcation surface (p — tp± which approach the observer with the speed of light at the retarded 
time. This is a non-integrable singularity. To superpose the fields of the constituent volume elements of an extended 
source we need to integrate the expression that appears on the right-hand side of Eq. (2) over the volume of the 
(r, 0, z) space occupied by that source. The multiple integral that needs to be evaluated here entails an integral 
with respect to (p whose integrand is proportional to \tp± — <p\~i in the vicinity of tp = tp± and so does not exist. 
Non-integrable singularities of this type commonly arise in the solutions of the wave equation over odd-dimensional 
space-times 0; they can be handled by the following method from the theory of generalized functions, a method 
originally devised by Hadamard [l6l ITtI Il8| . 

C. Hadamard's regularization technique 

Hadamard's regularization technique is applicable to situations in which the superposition of the potentials of the 
volume elements of an extended source yields a differentiable function of the observer's space-time coordinates, while 
the superposition of the fields of those same elements results in a divergent integral. This technique enables one to 
extract the physically relevant, finite value of the field of the extended source in question (which would follow from 
the differentiation of the potential owing to its entire volume) directly from the divergent integral that describes the 
superposition of the singular fields of its constituent elements. The calculation and subsequent differentiation of the 
potential owing to the entire source, a task which can hardly ever be performed analytically in physically realistic 
situations of this kind, is thus rendered unnecessary. 

In the present case, the singularity of the Lienard-Wiechert potential is like |1 — n • x/c| _1 ~ \(p± — <£>|~ 5 and so is 
integrable. The difficulty referred to above would not arise if we superposed the potentials, instead of the fields, of 
the consitituent volume elements of the source. We would know how to evaluate the required integral over (r, tp, z), at 
least in principle, and the function of (xp, tp) that we would thus obtain for the retarded potential of the entire source 



p 



FIG. 3: The bifurcation surface (i.e. the locus of source points that approach the observer along the radiation direction with 
the speed of light at the retarded time) associated with the observation point P at the observation time ip (the motion of the 
source is clockwise). The cusp Cb, along which the two sheets of the bifurcation surface meet, touches and is tangent to the 
light cylinder (f = 1) at a point on the plane passing through P normal to the rotation axis. This cusp curve is the locus 
of source points which approach the observer not only with the speed of light (dR/dt = — c) but also with zero acceleration 
(d 2 R/dt 2 = 0) along the radiation direction. For an observation point in the radiation zone, the spiralling surface that issues 
from P undergoes a large number of turns, in which its two sheets intersect one another, before reaching the light cylinder. 

could then be differentiated to obtain a finite, singularity-free expression for the field. The non-integrable singularity 
we have encountered stems from interchanging the orders of integration and differentiation, from differentiating the 
Lienard-Wiechert potentials of the individual source elements (to obtain the fields) prior to integrating them over the 
source volume. Though feasible in principle, it is not of course practical to calculate the retarded potential of any 
physically viable extended source of the type we are considering explicitly. The regularization method we are about to 
describe is such, however, that the finite value it would assign to the divergent integral we have encountered exactly 
equals the value of the field which would follow from directly differentiating the retarded potential of the extended 
source [H,[l3- 

All the essential features of the mathematical problem that we face when atempting to integrate the Lienard- 
Wiechert field (2) with respect to the source coordinate ip are illustrated by the following simple example. Consider the 
function f(x, a) = (a — x)~2. If we first evaluate the integral 1(a) = f(x, a)dx of this function with respect to x and 
then differentiate the result (I — 2a^) with respect to a we obtain the well-defined quantity dl/da — aT^ . But if we 
first differentiate f(x, a) with respect to a and then attempt to integrate the resulting function df /da = — ^(a — x)~ 5 
over x, we obtain 

J(a) = [ dxdf/da=-(a-x)- 1 *\ x=a + a-$ ) (7) 
Jo 

a quantity which is divergent. The value which Hadamard's regularization technique picks out as the finite part of 
the divergent integral J (a) is aT^ = dl /da. 

When applied to the more general form dx F(x)(a — x)~z of the above divergent integral, in which F(x) is a 
regular function, Hadamard's procedure consists of performing an integration by parts, 

pa 

I dx(a- x)~iF(x) = 2(o - x)-%F(x)\ 

Jo 

-2<z-iF(0) - 2 J Q a dx (a - x)-idF(x)/dx, (8) 

and discarding the term which is divergent. The remaining finite part [consisting of the last two terms in Eq. (8)] 
is the value which Hadamard's regularization assigns to this integral; it is the value one would obtain if one first 
evaluated —2 J ° dx F(x)(a — x)~2 and then differentiated the result with respect to a. 




FIG. 4: Close up of a segment of the cusp curve appearing in Fig. 3. The figure shows the section < z — 'zp < 5 of the light 
cylinder (f = 1) and the two sheets <f>± of the bifurcation surface (the locus of source points that approach the observer along 
the radiation direction with the speed of light at the retarded time) in the vicinity of its cusp curve Cb (the locus of source 
points that approach the observer with the speed of light and zero acceleration) for an observer who is located at rp — 3, 
(pp = 0. The cusp curve Cb is symmetrical with respect to the plane z — zp passing through the observation point P. The 
value Gj m of the Green's function Gj inside the bifurcation surface diverges on the inner sides of the two sheets <f>+ and <j>-. 
The value Gj ou of Gj outside the bifurcation surface undergoes a jump across the strip bordering on the cusp curve onto which 
these two sheets coalescence (in the limit Rp 3> 1). 

The integral arising from the superposition of the Lienard-Wiechert fields of the constituent volume elements of 
a super luminally rotating charge distribution has an integrand which, in contrast to that of the integral in Eq. (8), 
is singular within the domain of integration (rather than on the boundary of this domain). Although such cases are 
seldom encountered in physics, their treatment using Hadamard's regularization is well established in the theory of 
hyperbolic partial differential equations [ltl. Il7j. [In the case of Cerenkov emission from an extended source, where 
there are no contributions from the source elements outside the bifurcation surface (the inverted Cerenkov cone issuing 
from the observation point), the corresponding singularity occurs on the boundary of the domain of integration.] Here, 
we have a divergent integral whose Hadamard finite part consists of two terms, an integrated term which turns out 
to decay like Rp~^ with the distance Rp from the source as Rp tends to infinity, and an integral identical to the 
classical expression for the retarded field of a volume source [Eq. (14) of [Tjj] which decays spherically, like Rp^ 1 . 

D. The nonsphericallly decaying part of the emission 

The nonspherically decaying component of the radiation is detectable only within a limited region of space and 
during a limited interval of time: only when the cusps of the envelopes of wave fronts that emanate from the 
superluminally moving volume elements of the source propagate past the observer (see Fig. 2). 

The radiated field entails, at any instant during this limited time interval, a set of wave envelopes (each associated 
with the wave fronts emanating from a specific member of a corresponding set of source elements) whose cusps pass 
through the position of the privileged observer in question. These caustics arise from those volume elements of the 
source which approach the observer, along the radiation direction, with the speed of light and zero acceleration at 
the retarded time (Figs. 3 and 4). It is the contribution toward the intensity of the field from the filamentary locus 
of such source elements that decays like i?p _1 instead of Rp~ 2 . 

The nonspherical decay of the field does not contravene conservation of energy. The focused wave packets that 
embody the nonspherically decaying pulses are constantly dispersed and reconstructed out of other waves, so that 
the constructive interference of their constituent waves takes place within different solid angles on spheres of different 
radii Rp (see Appendix D of The integral of the flux of energy across a large sphere centred on the source is 
the same as the integral of the flux of energy across any other sphere that encloses the source. The strong fields that 
occur in focal regions are compensated by weaker fields elsewhere, so that the distribution of the flux of energy across 
such spheres is highly non-uniform and i?p-dependent. 
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Symbol 


Definition 




Angular rotation frequency of the distibution pattern of the source 




The angular frequency with which the source oscillates (in addition to moving) 


f 


Frequency of the radiation generated by the source 


n = 2nf/uj 


The harmonic number associated with the radiation frequency 


m 


The number of cycles of the sinusoidal wave train representing the azimuthal 
dependence of the rotating source distribution [see Eq. (9)] around the 
circumference of a circle centred on, and normal to, the rotation axis 


| Q ± mui | 


The two frequencies at which the nonspherically decaying component 
of the radiation from the source described in Eq. (9) is emitted. (The spectrum of the 
corresponding spherically decaying component of the radiation is limited to these two 
frequencies only if Q/ui is an integer.) 



TABLE I: Definition of the various frequencies and numbers used to describe the source and the emitted radiation 

E. Polarization of the source 

This paper is specifically concerned with the spectral and polarization properties of the nonspherically decaying 
component of the radiation from a superluminal source. To assist in identifying the origins of the various polarization 
components in the emitted radiation, we base the analysis on a representative polarization current j = dP/dt for 
which 

P r , v ,z(r, tp, z, t) = s rtVtZ {r, z) cos(mtp) cos(fit), -n < tp < n, (9a) 

with 

tp = tp — cut, (96) 

where P r ,ip,z are the cylindrical components of the polarization (the electric dipole moment per unit volume), s(r, z) 
is an arbitrary vector that vanishes outside a finite region of the (r, z) space and m is a positive integer. Equation (9) 
generalizes the earlier calculation reported in which was concerned with a rotating charge distribution, to a case 
in which the emitting polarization current flows in the r and z as well as in the tp direction and has a distribution 
pattern that oscillates in addition to moving. This is similar to the polarization currents available within recently 
constructed machines built to study the physics of superluminal emission (see H Elll E B0 and Appendix A 
of E). 

Note that the ranges of values of both tp and t in Eq. (9) are infinite, as in the case of a rotating point source, but 
the Lagrangian coordinate tp, which labels each source element by its azimuthal position at t = 0, lies in an interval of 
length 2-7T For a fixed value of t, the azimuthal dependence of the density (9) along each circle of radius r within 
the source is the same as that of a sinusoidal wave train, with the wavelength 27rr/m, whose m cycles fit around 
the circumference of the circle smoothly. As time elapses, this wave train both propagates around each circle with 
the velocity roj and oscillates in its amplitude with the frequency f2 (Table I). The vector s is here left arbitrary in 
order that we may investigate the polarization of the resulting radiation for all possible directions of the emitting 
current (Table II). Note that one can construct any distribution with a uniformly rotating pattern, P r . ¥ ,. z (r, tp, z), by 
the superposition over m of terms of the form s ri¥ , i2 (r, z, m) cos(mtp). 

III. DETAILED MATHEMATICAL TREATMENT OF THE PROBLEM 

A. Integral representation of the Green's function 

In the absence of boundaries, the retarded potential arising from any localized distribution of charges and 
currents with a density is given by 

A>*{xp,tp) = c- 1 Jd 3 xdtj tl (ii,t)5{tp-t-R/c)/R, fi = 0,---,3 (10) 

where 5 is the Dirac delta function, R stands for the magnitude of R = xp — x, and // = 1,2,3 designate the 
spatial components, A and j, of A 11 and j M in a Cartesian coordinate system. For the purposes of calculating the 
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electromagnetic fields 



E = - V P A" - dA/d(ct P ) and B = V P x A 



(11) 



generated by the source in Eq. (9), the space-time of source points may be marked either with (x, t) = (r, tp, z, t) or 
with the coordinates (r, tp, z, t) that naturally appear in the description of that rotating source. In fact, once tp is 
adopted as the coordinate that ranges over (— n, it), the retarded position ip of the rotating source point (r, tp, z) as 
well as the retarded time t could be used as the coordinate whose range is unlimited. 

The electric current density j = dP jdt that arises from the polarization distribution (9) is given, in terms of <p and 
tp, by the real part of 



J = l iw Eu=u+ /uexp[-i(ju<£ - n<p/w)]a, -it < tp < it, 



(12) 



where /i± = (Q/u>) ± m. Changing the variables of integration in Eq. (10) from (x, t) = (r, tp, z, t) to (r, tp, z, tp) and 
replacing j by the expression in Eq. (12), we obtain 



A = \i{w/c) Iv dV » exp(-w)s f A(p dp exp(iQ,tp/(j)8(g - 4>)/R(<p), 

where dV = rdrdtpdz. Here <fi stands for tp — tp P with tp P ~ tp P — Lutp, R(tp) is 

R(tp) = [(z P - zf + r P 2 + r 2 - 2r P r cos(p P - tp)] 3 , 

the function g is defined by 

g = tp-p P +R(tp), 



(13) 
(14) 
(15) 



with R = Ruj/c, Aip is the interval of azimuthal angle traversed by the source, and V is the volume occupied by the 
source in the (r, tp, z) space. 

Terms of the order of R~ 2 in VpX A, which do not contribute toward the flux of energy at infinity, may be discarded, 
as usual |l5j . if we are concerned only with the radiation field. Since the problem we will be considering entails the 
formation of caustics, however, we need to treat the arguments of the delta function and its derivative in the resulting 
expression for B exactly. Approximating \7p{R~ 1 5(t — tp + R/c)] by R~ 1 5'(t — tp + R/c)Vp(R/c), i.e. discarding the 
term that arises from the differentiation of R~ , we can therefore write the magnetic field of the radiation as 

B ~ ii(w/c) 2 E M=P± Iv dV f 1 ex P(" W) Ia v d( P exp(if^/w)n x s 

x5'(g-4>)/R{p), (16) 

where n = R(tp)/R(tp) and 6' denotes the derivative of the delta function with respect to its argument. 

To put the source density s = s r e r + s v e v + s z e z into a form suitable for inserting in Eq. (16), we need to express the 
(^-dependent base vectors (e r ,e v ,e z ) associated with the source point (r,tp,z) in terms of the constant base vectors 
(e rp ,e VP ,e 2p ) at the observation point (rp, tp P , zp): 



cos(p — tpp) sin(<yj — tpp) 
— sm(tp — tpp) cos(p — tpp) 
1 













- 6 2p _ 



Equation (17) together with the far-field value of n, 



lim n = sin#pe r „ 

P^oo 



cos9pe Zp , 9p — arctan(rp /zp), 
yields the following expression for the source term in Eq. (16): 



(17) 



(18) 



[s r cos dp cos(tp — tp P ) — s v cos dp sin(</5 — <pp) — s z sin 0p]&\\ 

+[s v cos(p - tp P ) + s r sin(p - <pp)]e±, (19) 



where = e VP (which is parallel to the plane of rotation) and ej_ = nxeii comprise a pair of unit vectors normal to 
the radiation direction n. 
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Inserting Eq. (19) in Eq. (16), rewriting S'(g — 4>) as —dS(g — (j))/dtp and making use of the fact that n, s and Ap 
are independent of tp, we arrive at 



B 



±i(w/c) 2 E M =^ ± Iv dVA*eoq)(-w) Ei=i UjdGj/dtp 



where 



ui = s r cos open + s v e±, U2 = — s v cosupen + s r ex; U3 = — s z smflpeii 



and Gj (j = 1, 2, 3) are the functions resulting from the remaining integration with respect to tp: 



G 2 
G 3 



ay — exp(il2y/w) 

it 



cos(y — (pp) 
sin(y — yp) 
1 



(20) 



(21) 



(22) 



The corresponding expression for the electric field is given by E ~ nxB, as in any other radiation problem. 

The functions Gi(r, tp, z; rp, pp, zp, pp) here act as Green's functions: they describe the fields of uniformly rotating 
point sources with fixed (Lagrangian) coordinates (r, tp, z) whose strengths sinusoidally vary with time. In the special 
case in which Q = 0, i.e. the strength of the source is constant, G3 reduces to the Green's function called Go in p| 
and represents the Lienard-Wiechert potential of the point source described in Eq. (1). This may be seen by noting 
that the evaluation of the delta function in Eq. (22) yields 



G, 



0=0 



E 

v=<pj 



1 



R\dg/d<p\' 



(23) 



where tpj are the solutions of the transcendental equation g{tp) = cj), solutions that are related to those of the equation 
h(t) = t + R(t)/c = tp for the retarded times via ip = (p + tot. The curves representing g(ip) versus (p have precisely 
the same forms as those appearing in Fig. 1. 

Similarly, the vector d(G2 cos8pe± + Gie^)/d(p is proportional to the Lienard-Wiechert field (2) when O is zero: the 
electric current associated with the rotating point source from which field (2) arises flows in the azimuthal direction 
and so corresponds to s r = s z = 0. The singularity structures of Gi are determined by the zeros of dg/dip or h [see 
Eq. (4)] and so are identical to the singularity structure already outlined in connection with Go in 



B. Asymptotic expansion of the Green's function in the time domain 

The retarded times at which the value of the Green's function (23) [or that of the Lienard-Wiechert field (2)] 
receives divergent contributions from the point source (r, tp, z) are given by the following solutions of dg /dp = [or, 
equivalently, of Eq. (4)] : 

<P± — VP + 27r - arccos[(l =p A^)/(ffp)], (24) 

with 

A = (f 2 P -l)(f 2 -l)-(z-zp) 2 (25) 

in which (f, Z] rp, Zp) stand for (ru)/c, zw/c; rpto / c, zpto/c). Note that for a given observation point (rp, (pp, Zp, tp), 
the critical times t± = (<p± — <p)/u> exist (i.e. are real) only when the (r, z) coordinates of the source point lie within 
the region A > of the (r, z) space shown in Fig. 5. 

The locus of source points in the (r, tp, z) space for which dg/dp — can be found by inserting t± = (p± — tp) /u) in 
the equation tp — t — R/c = which specifies the retarded times, or equivalently by inserting Eq. (24) in g — 4> = 
[see Eqs. (10), (13) and (15)]. The result is 

4, = 0± = g(p±) = 2tt - arccos[(l T A^)/(ff P )] + R±, (26) 

where 

R± = R(p±) - P - zp) 2 +r 2 +r 2 P - 2(1 T A*)]*. (27) 

The two sheets (±) of the tube- like spiralling surface tp — tp P 4- 4>±{r,z) meet tangentially and form a cusp where 
A = (Figs. 3 and 4). This cusp curve on which d 2 g/dp 2 (and h) as well as dg/dtp (and h) are zero, constitutes the 
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Z — Zp 



o 




r 



FIG. 5: The projection A = of the cusp curve of the bifurcation surface onto (r,z) space. The two sheets cj>±(r,z) of the 
bifurcation surface exist only for the values of (r, z) in A > 0. For a given (r, z) in A > 0, the source point (r, 0, z) lies within 
the bifurcation surface for (/>- < — &p < <t>+ an d outside this surface for other values of <p. The source elements whose (r, z) 
coordinates fall in A < approach the observer with a speed dR/dt < c at the retarded time and so make contributions toward 
the field that are no different from those made in the subluminal regime. The projection of the cusp curve of the envelope 
of wave fronts onto the (rp,zp) space is also given by A = and has the same shape: the function A is invariant under the 
transformation (r, z; rp, zp) — ► (rp, zp\ r, z). That the cusp curve of the envelope approaches the cone 9p = arcsin(l/f ) as Rp 
tends to infinity can be seen here from the slopes of the asymptotes of the curve A = 0: these asymptotes, and so the segments 
z > zp and z < zp of the cusp curve, lie in the plane z — zp for r = 1+ but open up and tend towards the vertical as f 
becomes increasingly greater than 1. 

locus of source points which approach the observer, not only with the speed of light, but also with zero acceleration 
along the radiation direction. On it, the coalescence of two neighbouring stationary points of the phase function 
[curve (a) of Fig. 1] results in a point of inflection [curve (b) of Fig. 1] and so in a stronger singularity of the Green's 
function. 

For source points in the vicinity of this cusp curve, a unifrom asymptotic approximation to the value of the integral 
(22) defining the Gi can be found by the method of Chester, Friedman and Ursell [20l l2ll | . Where it is analytic (i.e. 
for all x ^ xp), the function g(ip) in the argument of the delta function in Eq. (22) may be transformed into the 
following cubic function: 



where v is a new variable of integration replacing ip, and the coefficients c\ and C2 are chosen such that the values of 
the two functions on opposite sides of Eq. (28) coincide at their cxtrcma: 



g((p) = - ci 2 u + c 2 , 



(28) 



d = (f)*fa+-^-)*i C 2 = !(</>+ + </>_) 



(29) 



Insertion of Eq. (28) in Eq. (22) results in 




(30) 



where 




and Ais is the image of Aip under transformation (28). 

The leading term in the asymptotic expansion of the integral (30) for small c\ can now be obtained [2(], l22L I23L |24| 
by replacing its integrand fj with pj + qjp and extending its range to (— oo, oo): 




(32) 
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where 

Pj = |(/j'U=ci + fj\v=-cx), (33a) 



1j = 5 c i 1 (fj\v=a - fj\u=-n), (336) 

and the symbol ~ denotes asymptotic approximation. (Note that the critical points (p — ip± transform into v = =pci , 
respectively.) This integral has precisely the same form as that evaluated in Appendix A of 0. 

It behaves differently inside (0_ < (f> < (f> + ) and outside ((f> < <j>_, > <p + ) the bifurcation surface so that 

^-{g^ \ x \>i, (34a) 



with 



Gj m ~ 2ci 2 (1 — X 2 ) 2 [pj cos(| arcsinx) — ciOj sin(| arcsinx)], 



(346) 



Gj OUt ~ ci 2 (x 2 - 1) 2 sinh(iarccosh|x|) + ciq :) sgn(x) sinh(|arccosh|x|)], (34c) 

and x = 3(0 — C2)/(2ci 3 ) [cf. Eqs. (A16) and (A17) of 5]]. [The two-dimensional loci x = ±1 across which each Gj 
changes form correspond, according to Eq. (29), to the two sheets <fi± of the bifurcation surface, respectively] Explicit 
expressions for the coefficients Pj(r, z) and qj(r, z) are given in Appendix A. 

The above results show that as a source point (r, ip, z) in the vicinity of the cusp curve approaches the bifurcation 
surface from inside, i.e. as x —> 1 — or X —* — 1+j Gj m an d hence Gj diverges. However, as a source point approaches 
one of the sheets of the bifurcation surface from outside, Gj tends to a finite limit: 

Gr\^ ± = G/ Ut | x=±1 ~ ± 2c 1 a J )/(3c 1 2 ), (35) 

for the numerator of Gj OUt is also zero when \x\ = 1. The Green's function Gj is singular, in other words, only on the 
inner side of the bifurcation surface (see Fig. 6). 

Moreover, like the diffraction field near a focal point [22^ . the Green's function Gj OUt undergoes a phase shift 
across the coalescent surfaces cf> — <j>± at the cusp curve (Fig. 4 and 6). The shift in the sign of the second term 
in Eq. (34c) results in a finite discontinuity in the value of Gj OUt across the strip bordering on the cusp curve 
where the two sheets of the bifurcation surface are tangential: even in the limit c\ — ► 0, where <p- and 0+ are 
coincident, Gj OUt \ < j >= ^ )+ — G° nt \ ( j )=l j > _ has the non-vanishing value ^qj/ci. It is this discontinuity in the value of the 
Green's function (the potential of a point source) across the cusp curve of the bifurcation surface that gives rise to 
nonspherically decaying boundary contributions to the field of a volume source. (Note that a discontinuity of this 
type, which resembles that of a step function, cannot be handled by means of an analysis in the frequency domain 
unless the contributions from an infinitely large number of Fourier components of the decomposed Green's function 
are accurately superposed again when evaluating the field.) 



C. Hadamard's finite part of the integral representing the radiation field 

The integral which represents the superposition of the contributions of the source elements in the vicinity of the 
cusp curve of the bifurcation surface to the value of the magnetic field, i.e. the integral resulting from the insertion 
of Eq. (34) in Eq. (20), is divergent: the derivative dGj/dcp of the Green's function has a singularity [~ (1 — x 2 )~ 2 ] 
that is not integrable (with respect to x or equivalently (p). 

Our task in this section is to extract the Hadamard finite part p^.ll7| of this divergent integral, the finite quantity 
that we would have obtained had we been able to evaluate the potential A(xp, tp) of the entire source explicitly prior 
to applying the operator Vp x to this potential. 

The function dGj/d<p appearing under the integral sign in Eq. (20) is given by different expressions in different 
regions of (r, ip, z) space [see Eq. (34)]. For an observation point, the cusp curve of whose bifurcation surface intersects 
the source distribution, therefore, the domain of integration in Eq. (20) has to be divided into a part V\ n that lies 
within and a part V ov ± that lies without the bifurcation surface before the integrand of the integral in this equation 
can be written out explicitly (Fig. 3). The Green's function Gj has no singularities in the region A < of the (r, z) 
space (Fig. 5), so that the contribution of the source elements in the part of V ou t for which A < are no different 
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FIG. 6: The \ dependence of the Green's function Gi in the far-field limit where pi = 0. The values of this function inside and 
outside the interval —1 < \ < 1 represent Gi m and Gi out , respectively: the space inside the bifurcation surface [the interval 
4>- < <f> < (f)+ at a given (r, z)\ is mapped into \x\ < 1 and that outside the bifurcation surface into \\\ > 1- The function 
Gi m diverges on the inner sides, x = 1~ an d X — ~ l+i °f the two sheets (j>+ and <j>- of the bifurcation surface. The function 
Gi° ut is discontinuous across a two-dimensional strip bordering on the cusp curve of the bifurcation surface: even at points 
close to this cusp curve where the separation <f>+ — <f>- of the two sheets of the bifurcation surface tends to zero (Figs. 3 and 
4), the vanishingly small interval 0_ < <j> < <j>+ in <f> is mapped into the finite interval — 1 < x < 1 m Xj so that the difference 
Gi° ut | x =i+ — Gi out | x= _i_ between the values of the function Gi out on opposite sides of the strip in question remains non-zero. 

from the contributions of the source elements that lie in the subluminally moving portion of the source [see Eqs. 
(23)-(25)]. It would be sufficient to consider the contributions only of those source elements for which A > 0. All 
source elements for which dGj/dip is singular are taken into account once the integration with respect to tp (at fixed 
values of r and z in A > 0) is performed over the following two intervals: the interval 0_ < (p — ipp < <f>+ in which 
x| < 1 and Gj — Gj ln , and the remaining part of — n < <p < it in which \x\ > 1 and Gj = Gj OUt (see Fig. 6). 



Thus the magnetic field Ba>o that arises from the source elements in A > can be written, according to Eq. (20), 
as B in + B out with 



B 



in, out 



ii(w/c) 2 ^ 3 =1 J A > Q r dr dz UjKf 



in, out 



(36a) 



where 




(366) 



and 




(36c) 



(Recall that dV = rdrdipdz and <f> = (p — <fp.) 

Once it is integrated by parts, the integral in Eq. (36b) in turn splits into two terms: 




(37) 



of which the first (integrated) term is divergent [see Eq. (34) and Fig. 6]. Hadamard's finite part of K™ and hence of 
B ln (here designated by the prefix T) is obtained by discarding this divergent contribution toward the value of Kj m 
(see [HE El): 




(38) 
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Note that the singularity of the kernel of this integral, i.e. the singularity of Gj ln , is like that of \ip± — ip\~ s and so is 
integrable. 

The boundary contributions that result from the integration of the right-hand side of Eq. (36c) by parts are well- 
defined automatically: 

r r 1 ^+ 

K J OUt = T, M=li± { [^exp(-w)G/ ut J 

+i(/t-^+/;;^)#-" 2 ^p(-w)^° ut } 1 (39) 

for Gj° nt tends to a finite limit as the bifurcation surface is approached from outside (Fig. 6) and Gj |^=7r-y P equals 
Gj OUt |0 = _ 7r _ l pp when <f>± ^ ±n — (pp. The integral representing Kj OUt , in other words, is finite by itself and needs no 
regularization. 

If we now insert T{Kj in } and K° nt from Eqs. (38) and (39) in Eq. (36a) and combine B in and B out , we arrive at 
an expression for the Hadamard finite part of Ba>o which entails both a volume and a surface integral: J-{Ba>o} = 
B s + B ns . The volume integral 

B s = \i^lcf E M=M± f J A >o rdr dz # exp(- w) E?=i u j^j (40) 

has the same form as the familiar integral representation of the field of a subluminal source an d decays spherically 
(like Rp^ 1 for Rp — > 00). The surface integral 

B- = -hKu/crZUfA^rdrdzuiK^™^ (41) 



stems from the boundary contribution 



^■ boundary _ 



Uxp(-w)G/ ut f + (42) 

L J ch_ 



in Eq. (39) and, as we shall see below, turns out to decay nonspherically (like Rp 2 for Rp — > 00). 
Making use of Eq. (35) to rewrite Gj OUt |^+ in Eq. (42) explicitly, we obtain 



^.boundary = 2 Me x P [-i M (^ + C 2 )] [2^"^ COs(fW) 

-ici-^sind^ci 3 )], (43) 



where ci and C2 arc defined in Eq. (29). The asymptotic expansion of Gj OUt in Eq. (34c) is for small ci. To be 
consistent, therefore, we must likewise replace the above expression by the first term of its Taylor expansion in powers 
of c\: 

^.boundary = 4^-1^. ^[-^P + <f>-)} +'", (44) 

for the remainder of this series is by a factor of the order of ci 2 smaller than the above retained term. 
The value of c\ close to the cusp curve of the bifurcation surface (where A = 0) is in turn given by 

c 1 = 2-^(f 2 f 2 ;-l)-^A5+0(A). (45) 

This may be obtained by using Eq. (25) to express z everywhere in Eqs. (26) and (27) in terms of A and f and 
expanding the resulting expressions in powers of As. When the observation point lies in the far zone, f on the cusp 
curve of the bifurcation surface has the value csc#p [see Eq. (25)] and so c\ assumes the value 2 _ 3i?p 1 A2 . Moreover, 
according to Eqs. (A9) and (A10), 

qj ^2^(uj/c)R p 1 exp[i(n/uj){(p P + 3n/2)}q 3 (46a) 

with 

g-j = (1 - ifi/w iO/w) (466) 



for A <C 1 and Rp 3> 1 (see Appendix A). 
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Equations (41), (44), (45) and (46) thus jointly yield 

B ns ~ -|i exp[i(n/w) (cpp + 3rr/2)] £„ =/i± M exp(- Wp) 

x Ei=i % I A >o fdfd2 A ^ u i exp(-i/i0_). (47) 

The integrand of the surface integral in Eq. (47) is still singular on the projection of the cusp curve of the bifurcation 
surface onto the (r, z) space (Fig. 5) but this singularity is integrable. 

D. Amplitude, frequencies and polarization of the nonspherically decaying radiation outside the plane of 

source's orbit 

The leading term in the asymptotic expansion of the integral in Eq. (47) for high ji± (i.e. for the regime in which 
the wavelengths of oscillations of the source are much shorter than the length scale c/u> of its orbit) may be obtained 
by the method of stationary phase provided that Op ^ w/2. 

Both derivatives, d<f>-/dr and d(j)-/dz, of the function that appears in the phase of the rapidly oscillating expo- 
nential in Eq. (47) vanish at the point f — 1, z = Zp, where the cusp curve of the bifurcation surface touches, and 
is tangential to, the light cylinder (see Figs. 3 and 4). However, d 2 <p-/dr 2 diverges at this point, so that neither the 
phase nor the amplitude of the integrand of the integral in Eq. (47) are analytic at f = 1, z — zp. Only for an observer 
who is located outside the plane of rotation, i.e. whose coordinate zp does not match the coordinate z of any source 
element, is the function 0_ analytic throughout the domain of integration. To take advantage of the simplifications 
offered by the analyticity of 0_ as a function of f, we proceed under the assumption that Op =^ tt/2. 

Since 0_ = g{ip~) and dg/d(p — dg/dr — along the curve 

C : r = f c (z) = {\{rl + 1) - [\(f P - l) 2 - (z - z P ) 2 ]*}i 

tp = (fc(z) = <pp + 2-7T — arccos(fc/rp) (48) 

[Eqs. (15) and (26)], 0_ is stationary as a function of f, i.e. d(f)-/dr = 0, on the projection fc(z) of C onto the (r, z) 
plane (see also [lj|)- In the far-field limit, where the terms [z — zp) 2 /(fp — l) 2 and rc/rp in Eq. (48) are much 
smaller than unity, curve C coincides with the locus 

C b : r=[l + (z-z P ) 2 /(r p -l)]K 

ip = (pp + 2tt — arccos[l/(ffp)], (49) 



of source points which approach the observer along the radiation direction with the wave speed and zero acceleration 
at the retarded time, i.e. coincides with the cusp curve of the bifurcation surface (Figs. 3 and 4). 

It can be seen from the far-field limit of Eq. (49) that the cusp curve Cb would intersect the source distribution 
if Op lies in the interval | Op — -| < arccos (l/f>), where r> = r>cj/c and r> is the radial coordinate of the outer 
boundary of the source [5j. Hence, the range of values of Op for which the stationary points of 0_ fall within the 
source distribution would be as wide as (0, w) if the extent r> — c/uj of the supcrlumially moving part of the source 
is comparable to the radius c/u of the light cylinder (see Fig. 5). 

The dominant terms in the Taylor expansion of 0_ about a point (fc, <f>Ci z) on curve C (with an arbitrary 
coordinate z =/= zp) are 

4>-=4>c- \Wl - IWI - I)" 1 - l]Rc\f - fc? + ■■■, (50) 

in which we have denoted the values of <fi and R on C by <pc = Rc + fc ~ <Pp an d 

R c = [(zp - z) 2 + fp - f 2 c }i , (51) 

respectively. Moreover, has the finite value r 2 -, — 1 on C [Eqs. (25) and (48)]. 

We are now in a position to evaluate the leading term in the asymptotic expansion of the r integral in Eq. (47) 
by applying the principle of stationary phase j2j 13 : by replacing the phase of the rapidly oscillating exponential 
that appears in this integral with its Taylor expansion (50) and approximating the amplitude of the integrand with 
its value along C. The resulting integral, at a given value of z, is 

/ fdf A~ 2 Uj exp(— i/i0-) ~ exp(— in4>c)rc {Pc ~ l)~ lu j \c 

JA>0 
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x g>' fc dr, exp{i/x[i(f F ~ MR^Wh (52) 

in which fc, 0c and i?c have their far-field values 

fc—cac9p, (pc — 3ir/2 + Rc, Rc — R P — zcos8 P , (53) 

and T] = r — ?c [see Eqs. (48) and (51)]. This integral entails an integrand whose phase is large on account of both 
fi± 3> 1 and Rp 3> 1. It can be cast into the form of, and evaluated as, a Fresnel integral with a large argument to 
arrive at 

f i i - —i i 

/ fdf A~2 Uj , exp(— i[i(f)-) ~ (7r/2)2i? p 2 1 sin^p cos^pj -1 !^! - ^ 

Ja>o 

x exp[-i(/^ c - f sgn/i)]uj| c (54) 



when 6>p 7^ tt/2. Note that the values Uj|c of Uj along the curve C are functions of z and that this curve becomes 
parallel to the z axis as Rp tends to infinity (Figs. 4 and 5). 

Once the surface integration in Eq. (47) is expressed in terms of a double integral, this result may be used to obtain 

B ns ~ -|i(2^)5i? p ^| s in2^p|- 1 exp(ir!^ c / W )E AI=A1± M*sgn(/0 exp(if sgn/x) 

x exp[-i/i(_Rp - ustp + (p c )] J2j=i ?j X!^, ^ Uj| c exp(i^l cos 6 P ), (55) 



in which (pc = </?p + 3tt /2. The remaining 2 integration in this expression amounts to a Fourier decomposition of the 
source densities s r>(p>z \c- 

Insertion of Eqs. (21) and (46b) in Eq. (55) and the introduction of the Fourier transforms 

/oo 
dz s r ^, z | c exp(i^£ cos 6 P ). (56) 
-00 

result in 

E ns ~ |(2^)ii?;*|sin20p|- 1 exp(i^ c / W )E A1= , I± |/i|5 S gn(/x) exp(i| S gn/i) 
x exp[-i^i(ijp — utp + (pc)]{(istp + Qs r /w)e|| 
— [(is r — VLs^/uj) cosOp + Os z sin^p/cjje^} (57) 



for the electric field (E ns ~ nxB ns ) of the nonsphcrically decaying part of the radiation. This expression is valid, 
of course, only for an observation point, outside the plane of rotation, the cusp curve of whose bifurcation surface 
intersects the source distribution, i.e. for < \6 P — ?| < arccos (l/f>). 

Thus the spectrum of the nonspherically decaying part of the radiation only contains the frequencies h±<jJ = fiimw, 
i.e. the frequencies that enter the creation or practical implementation of the source described in Eq. (9). This contrasts 
with the spectrum of the spherically spreading part of the radiation which extends to frequencies of the order of f2 3 /w 2 
when Q/uj 3> 1 (cf. [Tfll]). The intensity of the radiation at the frequency Q + mu> is the same as its intensity at 
il — moj when either fi/w 3> m or Vl/u) <^ m. The radiation at one of these frequencies would have a much lower 
intensity, on the other hand, if f2 and muj are comparable in magnitude. 

Equation (57) shows, moreover, that the nonspherically decaying component of the radiation is linearly polarized 
both for Q/uj <C 1 and for il/u> 3> 1 when one of the cylindrical components of s is appreciably larger than the others. 
In the case of an s which lies in the azimuthal direction, the radiation is polarized parallel to the plane of rotation 
for CI/oj <C 1 and normal to that direction (and phase shifted by tt/2) for Q,/uj 3> 1. The plane of polarization of the 
radiation coincides with the plane passing through the observer and the rotation axis if s lies parallel to the rotation 
axis. For this radiation to be clliptically polarized, on the other hand, s needs to have two cylindrical components 
that are comparable in magnitude and 0,/lo has to be of the order of unity (Table [HJ). 



IV. DISCUSSION: COMPARISON OF THE NONSPHERICALLY AND SPHERICALLY DECAYING 

COMPONENTS OF THE RADIATION 



The radiation field that arises from the super luminal portion (r > c/uj) of the volume source described in Eq. (9) 
consists, as shown in the preceding section, of two components: a nonspherically decaying component whose intensity 
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% > 1 


n ^ 1 

UJ 


S<1 


Sr / 0, S v = S z = 


linear, en , phase= 


elliptic 


linear, ex, phase= — ^ 


Sip 0, S r = S z = 


linear, ex, phase= 


elliptic 


linear, ey, phase= +5 


S 3 / 0, S r = S v = 


linear, ex, phases n 


linear, ex, phases n 


linear, ex, phase= ir 



TABLE II: The state of polarization of the nonspherically decaying component of the emitted radiation for different ranges of 
Q/u) and different orientations of the emitting polarization current 

diminishes like Rp^ 1 with the distance Rp from the source and a spherically spreading component, one whose intensity 
has the conventional dependence Rp~ 2 on Rp. The former is described by Eq. (57) and the latter, which follows from 
Eq. (40), was earlier calculated in [l^. The nonspherically decaying part of the radiation E ns is only emitted at the 
two frequencies [a±u>, whereas the spherically decaying part E s has a discrete spectrum, comprising multiples nuj of 
the rotation frequency, which extends as far asn~ (f2/w) 3 when O/w (3> 1) is different from an integer (see Table I). 

As in the case of any other linear system, the present emission process generates an output only at those frequencies 
which are carried both by its input (the source) and its response (Green's) function. The agent responsible for the 
generation of a broadband radiation from the source described in Eq. (9), whose creation or practical implementation 
only entails the two frequencies /i±w, is acceleration: a remarkable effect of centripetal acceleration is to enrich the 
spectral content of a rotating volume source, for which Q/lo is different from an integer, by effectively endowing the 
distribution of its density with space-time discontinuities. (For a detailed discussion of this point, see [l^.'l 

The electric field of the spherically spreading part of the radiation for incommensurate values of £1 and uj is given 
by the real part of 

oc 

E s = E^ + 2^E>xp(-m^ P ), (58) 

n=l 

in which the Fourier component E^ of this field at the frequency nu has the following value beyond the Fresnel zone: 

E s „ ~ ifp 1 cxp{-i[n(i? P + |tt) - (Sl/u)(<pp + §7r)]}(f> - f<)Q^Q z 

+{m -> -m,Q -> -ft}, (59) 

with 

Q<p = - ^2 ^ 2 sin[7r(n - n)]/(n - fi), (60) 

and 

Q z = [s r J n -n/u>(n) +is v ,J' n _f2/ w (n)]e|| + [(s v cos9 P 

—s z sm6 P )J n _ n / u (n) -is r cos0pJ' n _Q /a; (n)]e_L (61) 

[cf. Eq. (66) of (l9|]. Here, f< < 1 and r> > 1 denote the lower and upper limits of the radial interval in which the 
source densities s ryVyZ are non-zero, the symbol {m — » — m, ft — * —ft} designates a term like the one preceding it in 
which m and ft are everywhere replaced by —m and —ft, respectively, and J n _n/w(?'0 &ud 3' n -Q/ui{ n ) are the Anger 
function and the derivative of the Anger function with respect to its argument. 

To compare the amplitudes of E ns and E s at a frequency with which both these components are emitted, let us 
consider a case in which fi/u> equals an integer, so that the spherically decaying part of the radiation is also emitted 
only at the frequencies n±u>. In this case, the quantity in Eq. (60) is non-zero only if n equals fi+ or and 

Qip\n— fi± 7^M± • 

At the higher of the two frequencies, i.e. at /i + oj = ft + raw, the amplitude of E ns has the value 

|E ns M+ | ~ l(2n)iR-* | sin20p|-V+5 1 (iS v + StS r /ui)e { \ 

— [(is r — ftstp/oj) cosOp + fts z sin^p/o;]ex| (62) 
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according to Eq. (57). This should be compared with the following amplitude implied by Eqs. (58)-(61): 

~ i 7 "V 1 (*'> -r<)n+ 2 \[s r J m {n+) +is v J' m (M+)]e|| 

+ [{s v cosOp — s z sin^p) J m (/i + ) — is r cos6*p J' m (/i + )]ej_| (63) 

since the Anger functions J n -n/uj( n ) an d J' '„_n / 'iu( n ) m Eq. (61) respectively reduce to the Bessel functions J m .(/i+) 
and J' m (/i+) when Cl/ui is an interger and so n is exactly equal to /i+ (see jl9i l25j V 

In a case where the emitting polarization current is parallel to the rotation axis, for instance, the ratio of the 
amplitudes of the two components of the radiation is given by 

lE-VI/lE^l^^^ilsec^K-i^MlJ™^)!- 1 ^^^-^)- 1 , (64) 

since s r and s v would then be zero. For fi/a> ^ 1, the amplitude of J m (M+) is of the order of Irrespective 
of whether m is smaller than or comparable to Q/uj, therefore, the factor (Q/uj)\ J m (/i+)| _1 is independent 

of frequency. Thus the above ratio is already much greater than unity (~ /i + 5 ) at the Fresnel distance Rp ~ 
(/i + cj/c)(r> — r<) 2 from the source, a result which holds true, as can be seen from Eqs. (57) and (59), even when 
Q/lu is different from an integer and s r and s v are non-zero. 

V. SUMMARY 

We have examined the electromagnetic emission from polarization charge-currents whose distribution patterns have 
the time dependence of a travelling wave with an accelerated superluminal motion. Such macroscopic polarization 
currents are not incompatible with the requirements of special relativity because their superluminally moving distri- 
bution patterns are created by the coordinated motion of aggregates of subluminally moving particles. Our analysis 
is based on an emitting polarization current, which has a poloidal as well as a toroidal component, and which has a 
distribution pattern that sinusoidally oscillates in addition to superluminally rotating (Table I); similar currents are 
employed in recently-constructed machines designed to test the physics of superluminal emission. 

We find that such sources generate localized electromagnetic waves that do not decay spherically, i.e. that do not 
have an intensity diminishing like Rp 2 with the distance Rp from their source [Eq. (57)]. The nonspherical decay of 
the focused wave packets that are emitted does not contravene conservation of energy (Section IID): the constructive 
interference of the constituent waves of such propagating caustics takes place within different solid angles on spheres 
of different radii (Rp) centred on the source. 

Detailed analysis in the far-field limit shows that the spectrum of the nonspherically decaying part of the radiation 
emitted by the source described in Eq. (9) only contains the frequencies Q. ± mu, i.e. the frequencies that enter the 
creation or practical implementation of that source (Section HID and Table I). This contrasts with the spectrum of 
the spherically spreading part of the radiation which extends to higher frequencies (Section IV) . 

We have also determined the polarization of the nonspherically decaying component of the radiation in the far-field 
limit. In many cases, the emission is highly linearly polarized [Eq. (57) and Table II]; however, with certain source 
frequencies and polarizations it is possible to also produce elliptical polarized emission. 

Finally, we have examined the relative amplitudes of the spherically and nonspherically decaying components of 
the emission. We find that even at relatively short distances from the emitter, the latter component can represent 
the greater part of the observed signal. 
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Appendix: evaluation of the leading terms in the asymptotic expansion of the Green's function 

This appendix concerns the evaluation of the coefficients Pj(r, z) and qj(r, z) in the asymptotic expansion (34) of 
the Green's function Gj for small c\. These coefficients are defined, by Eqs. (31) and (33), in terms of the functions 
V±; R± an d ci, which appear in Eqs. (24)-(27) and (29), and the derivative dip/dis which is to be calculated from 
Eq. (28). 
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We have already seen in Section IIIC that close to the cusp curve of the bifurcation surface (where A = 0) the 
function C\(r, z) can be approximated as in Eq. (45). In the regime of validity of the asymptotic expansion (34), where 
A is much smaller than [f 2 r 2 p — 1)5 , the functions <p± and R± may likewise be expressed in terms of (r, A) [instead 
of (f, z)\ and expanded in powers of Aj to arrive at 

p>± = p cT (f 2 f 2 P - l)"i A* + <9(A), (Al) 



R± = (f 2 f 2 P - 1)4 ± (f 2 f| - 1)-*A* + 0(A), (A2) 

where p c = p P + 2ir — arccos[l/(rf P )]. Recall that v — ±Ci in Eq. (33) are the images under the mapping (28) of 
p = p^:, respectively. 

The remaining functions dip/ dv\ v —± Cl that appear in the definitions of pj and qj are indeterminate. Their corre- 
sponding values have to be found by repeated differentiation of Eq. (28) with respect to v. 

(dg/dp)(dtp/dv) = v 2 - Cl 2 , (A3) 



(d 2 g/dp 2 )(dp/dv) 2 + (dg/dp)(d 2 p/dv 2 ) = 2v, (A4) 

etc., and the evaluation of the resulting relations at v = ±c\. This procedure, which amounts to applying the 
l'Hopital's rule, yields 

dp>/dv\ v=±ci = (2 Cl R T )?/Ai (A5) 

and 

d 2 y/dv 2 \ v =± Cl =±|A-i(2A F /ci)3[l-(2sciA F /Ai)i(l±3A3/^:)] {A6) 

for the values of the first two derivatives of tp. Close to the cusp curve A = 0, where c\ and R± may be approximated 
as in Eqs. (45) and (A2), these reduce to dpjdv\ v= Q = 2s and d 2 pjdv 2 \ v= Q = —2~*{r 2 r 2 P — 1)~2, so that the 
contribution of d 2 <p/ dv 2 \ u=0 to qj is negligible in the far field. 

Inserting Eqs. (Al), (A2), (A5) and (45) in the defining equations (31) and (33), keeping only the dominant terms 
in powers of A 2 and Rj. 1 , and taking the far- field limit Rp 3> 1 of the resulting expressions, we obtain 

Pi ~ 2^(uj/c)R~ 2 e7ip(iflp c /uj), (A7) 



p 2 ^ -23 (uj/cjkp 1 exp(irip c /u;), p 3 ~ -p 2 , 



(A8) 



and 



q\ ~ 23 (uj/cjRp 1 exp(iflp c /Lu), 



(A9) 



q 2 ~ -q 3 =± -\{Q,/uj)q u (AIQ) 

where ip c in these expressions has its far- field value pp + 3ir/2. Here, we have made use of the fact that the 
indeterminate ratio (A 2 /ci)a=o has the value 2 a (f 2 fp — 1)2 , and that as Rp tends to infinity f approaches the value 
cscOp along the cusp curve of the bifurcation surface [see Eq. (49)]. 
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